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Abstract 

A system of strongly interacting fermions in a solid state is discussed. 
As an example, a structure of singlet and triplet coupled 2-particle states 
and their excitation spectra are investigated. It is shown that an account 
of intersite fermion interaction leads to new boson modes. Their number 
depends on dimensionality of a system. It is shown also that a reasonable 
choice of Hamiltonian parameters leads to a small effective mass of local 
pairs. A problem of interaction of coupled states with external fields are 
discussed. The calculations are carried out for intermediate and strong 
fermion - fermion coupling. The same analysis can be performed for 
another sorts of pairings like d-wave. 



1 Introduction 



The discovery of high-temperature superconductors (HTSC) that differ drasti- 
cally from conventional superconductors leads to a great number of theoretical 
models. A small coherence length in HTSC-materials indicates that a size of 
fermion pairs, which is crucial for the superconducting properties, should be of 
the order of lattice parameter. 

In this paper the formalism of tightly binding fermions in solids is developed. 
The main tool is extended Hubbard model with a fermion-fermion intersite 
attraction. The origin of the attraction could be 'phononic', 'excitonic' Q, 
'plasmonic' [pi or 'magnonic' p|]. 
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Last years indicate the power of the bipolaron theory , |^ and one could 
identify the local pairs and bipolarons ~ local pairs with phonon mechanism of 
2-fermion attraction. The presence of the attraction changes the band structure 
significantly. Fermions being placed on different sites are coupled into pairs. 
The continuum of electron-hole pairs is situated above the local pair bands and 
does not overlap with them both in the normal and superconducting phases - 
see discussion in It is the main difference of BCS-like theories and local pair 
ones. 

The pair size seems to be of the same order of the distance between pairs 
that appears to be rather difficult for theoretical analysis. Nevertheless there 
are some papers Q, 0] that are devoted to this complex intermediate region 
of crossover from Cooper pair regime to local pair one. 

A presence of a superconducting state indicates an existence of boson type 
excitations, that are able to form a condensate at low temperatures. But there 
are fermions (electrons) in solids as the only charge carriers from the beginning 
that interact with phonons and other quasiparticles in crystal. This interaction 
in HTSC - materials is very complicated both for experimental and theoretical 
investigations. As a result, a large variety of superconductivity models has been 
developed. Most of papers deal with a ground state and excitation spectra on the 
basis of an attractive or repulsive Hubbard model (see review [^). Although 
a generalization to a more realistic and complicated situation is not obvious, 
several interesting papers have been recently appeared Q, In |^ 

the authors considered in detail a ground state of such a system with rather 
general form of interaction in site representation. In pO[ | the Green function 
method was applied to the boson spectra excitations. But the latter work was 
motivated by Auger experiments in noble metals and only fermion - fermion 
repulsion was taken into account. 

This work is devoted to the investigation of the boson excitations - their 
band structure, interaction with each other and with external fields. 

We start with rather general form of fermion - fermion interaction to be able 
to describe a large variety of interaction forms in different HTSC - materials. 
To be more definite, we postulate the effective Hamiltonian for fermi excitations 
in site representation in the form: 

H = cjjCM + ^ t{m - n)cliCN + 

M M,N 

T ^ UMN.M'N'C'^.fC'^CM'CN' + ''^ BmncIjCn- (l) 
MNM'N' MN 

Here ^ is the chemical potential, and for example M — (m, s) , where m rep- 
resents the site index, s is the spin projection of the fermion, c^^,cm are the 
operators for fermion creation and annihilation on the site m with the spin 
projection s = ±i, Umn,m'N' - is the matrix element of fermion interaction. 
This matrix is antisymmetric in its first and second pairs of indices separately: 
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Umn,M'N' = -Unmm'N' = ~Umn,N'M'- A hopping integral t(m - n) char- 
acterizes the fermion motion and describes a kinetic energy in a momentum 
representation. Bmn represents an arbitrary external field. 

It must be stressed here that only fermi statistics of primordial quasiparti- 
cles is important. One can keep in mind polarons, holes, magnons and so on. 
The question is if the interaction is strong enough to form the local pairs or 
we have Cooper pairs in high-Tj, superconductors as in the ordinary ones. A 
mathematical method in the latter case has been developed quite accurately. 
The method, developed in this paper for the description of the local pairs is 
based on a functional integral technique , that was widely used in different 
theoretical investigations (see, for example, It was applied 

to a single-polaron problem jl^ and is used here for the local pair description 
of strongly interacting fermions. Calculations appear to be rather simple in 
the case of short - distance interaction. Furthermore mathematical difficulties 
increase slightly with a number of accounted sites. 

Below the general formalism is presented for an arbitrary form of fermion - 
fermion interaction. Detailed calculations of spectra excitations and a structure 
of the pairs are discussed. Several applications of the developed method were 
considered in papers devoted to tunneling and where it was shown 
that experimental data on neutron scattering and infrared absorption spectra 
are in good agreement with predictions of local pair model. 



2 General formalism 
2.1 Main equations 

The Hamiltonian H in the form is quite general and its origin is not discussed 
here. The external field B is written in a general form too. This form includes 
the interaction with electromagnetic field and neutrons, for example. Let's start 
with partition function of the system in a functional integral form : 



Z — Spe 





A 




dT[^^C*ACA + ^^rABC*ACB - ^ jUAB-A'B'C*ACgCA'CB'], (2) 



A AB ABA'B' 



TAB = [P-- ^^pa)^AB- 
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Expression (|^) is written in momentum representation where complex index 
A — {pa, sa) includes quasimomentum and spin projection of a fermion; ujp^ is 
the fermion spectrum determined as a Fourier transform of the hopping integral; 
c\{t), ca{t) are independent grassman variables; /3 = 1/T is the inverse tem- 
perature. The external field B is omitted in this section. To get rid of the four 
grassman variables product in the exponent and introduce the bosonic variables 
V'AB (''") one can use the Hubbard - Stratonovich transformation ||l7|| , which 
is usually applied to consideration of bosonic fields of one argument that de- 
scribes boson motion as a whole. In this latter case all the information about the 
internal structure of the bound state is lost, though the results of conventional 
BCS - theory of superconductivity are reproduced But as far as we 

intend to study the internal pair structure and the interaction of the pairs with 
external fields we should use more complex bosonic fields V'ab(t), -(AyiB (■'')*■ In 
this case expression (H) is transformed into: 



Z = I Dc* DcD^j* Di^e"" , 



f 

/ dT[Y^C*ACA + Yl 
•^0 A AB 



rABC*ACB + -^i'AB^^ACB + -^i^ABC^Cs 



'4'ABUAh;A'B'^A'B']- (3) 



ABA'B' 



Our central problem is to obtain equations for the wave functions of the 
coupled states and to find their solutions for the chosen form of interaction 
term. An imaginary "time" Fourier - transformation gives: 



Se^Y^'^^B +i£SAB)cA{s)cB{e) + ^ ^ tp*As{e + e')cA{e)cB{e') 

eAB ee'AB 

+ 1 J2 i^ABie + e')c\ie)c*Bis') + J2 ^*ABiE)UXh-A,s,^A'B'{E), 



2 

ee'AB EABA'B' 



where Matsubara frequencies e — {2n + l)Tr/f3, e' = (2n' + l)7r//3 due to antiperi- 
odic boundary conditions for fermionic variables and E ~ 2Tm/ (3 due to periodic 
boundary conditions for boson dynamic variables. Now one can integrate out 
the grassman variables |12[: 



Z = J D^:* Dii^e^" , 

Sb^ E rAB{E)Ull^j,,s4A'B'{E)+Tr{lnY), (4) 

ABA'B'E 



4 



■y ( '\ ^ ( (tab +ie5AB)5e.e' iPab{£ + s') 

yAB{e,e)-y ^'^^(e + e') _^rAB+ie5AB)5,,e 

Boson dynamic variables i^^^jipAB having been introduced, the expression 
(H) remains not very suitable for the further analysis. Next step is to expand 
the action Sb around a classical action Sh[^ci\- {^Sb/Sipci = by definition). 
The simplest case ipd = 0, takes place: a) in a normal phase (by definition), b) 
near critical temperature, where a condensate density is small, and c) in a low 
fermion density limit. Calculations for ^jjd ^ were discussed in [^]. Here we 
will be interested in the normal phase of the system. 

Let's expand the effective action Sh (jj) in Tailor series up to the second 
order in fields V'ab'V'as around trivial classical trajectory t/jc; = 0. A fourth 
order in the fluctuation describes interaction of pairs and was discussed in po| . 
As it was shown in ref. this term is small in a low-density region and/or 

for a weak fermion interaction. Here we neglect the pair interaction and limit 
ourselves to the study of boson excitation spectra and their interaction with 
external fields. Simple, but tedious calculations lead to quadratic part of the 
effective boson action: 



Sf^ = Y.ABA'B'ErAB{E)bAB,A'B'{E)i,A'B'{E\ 

~ 1 fi(/3nA/2)+t/i(/3nB/2) 1 

Dab,A'B'[E) - J ■ —[E+TUTTl^"'^^'"^'^ + '^A'B';AB ' 

^A = ^ M- 

To remove the denominator in the first term and the inverse power of the 
interaction operator U, one has to make the following substitution of variables 
in the functional integral (Q): U~^ip i^^ipAB ' {^EQa - ^b) ^ab^ ■ This 
procedure leads to the partition function in the form: 

Z= D^*D'lljexp{ ^ i^AB{~ElAB,A'B' +iESAA'SBB')'fpA'B'}, (5) 
ABA'B'E 



where 



DaB,A'B' = i^A + ^b)SaA'SbB' + Wab,B'A', (6) 



WaB;B'A' = i ■ {th^ + th^)UAB:,B'A'. (7) 
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Now we have to diagonalize the quadratic form of the exponent, that is to solve 
the equation: 



DaB-A'B'U^A'B' = ^nU^AB- (8) 

A'B 



Complex index n numerating the eigenstates of the system will be shown to 
include quasimomentum of the local pair, its spin and a band number. The 
eigenvalues A„ obey the secular equation: 

Det{DAB;A'B' - KSaA'5bb') = 0. (9) 
The eigenfunctions are chosen to be orthonormal: 



E 

AB 



Snn'. (10) 



The final step consists of introducing new variables Cn{E) in the integral 

^AB = J2CniE)u'^Ah rAB = 



This leads to the expression: 

Z = l[ [ DClDCn exp{C:(^£; - A„)C„}. (11) 

This form of the partition function clarifies the meaning of the new variables 
C^,Cn as the wave functions of bosonic excitations of energy A,i. It will be 

(n) 

shown below that it^^ functions describing the internal structure of the bound 
state must be taken into account in the analysis of the pair interaction with 
an external field. Eq.(||) is the only one we need to calculate bosonic - type 
spectra excitations. The pair energy depends on its quasimomentum and spin 
projection. Quasimomentum and spin conservation can be accounted for in the 
form of interaction term: 

Wab.b'A' = \wAB.B'A'SiA + B-A'-B') (12) 

Index A — {pa,sa) includes both the momentum p^i and the spin projection 
SA of the fermion. Let variables R — A + B, R' — A' + B' characterize the 
momentum and spin projection of the coupled state as a whole. One can easily 
find from Eqs.([^,[T^), that matrix D is diagonal in indices R, R'. Several energy 
bands may exist and we specify them by additional index a. It means that the 
complex index n in the expression ( [ll| ) should incorporate the quasimomentum 
Pk of the 2-particle state as well as its spin projection = 0, ±1, united in 



6 



the common index K = (jjktCtk), and eigenfunctions of the matrix D should 
have the form: 

^AB — ^RA — ^RA — OK,RXjfj^ , 

where k = 0(1) corresponds to singlet (triplet) case. The substitution of ex- 
pression ( p^ into Eq. (||) along with the normalization condition term ( |lo|) lead 
to: 



[^A + ^K-A - K{K)]xk'a + \Y.^A'^-a-^k-A',A'^ka' = 0, (14) 

A' 

and 

KA KA — ^KK'^aa'- } 

A 

These equations determine the pair structure a^^^^ and the excitation spec- 
tra \K,,a{K) The previous analysis shows that for k = the solution is antisym- 
metrical in the fermion spin indices (x^'^^ = —x^k'k-a) 'v^hile it is symmetrical 
for K — 1. It confirms the notation k = for the singlet pairs and k = 1 for the 
triplet ones. 



2.2 Equations for triplet and singlet states 

Cooper pairs in the BCS theory are singlets. There are many reasons to think 
that in HTSC - materials triplet pairs exist along with singlet ones Both 
of these possibilities are included in the system of equations ( p^ ) and the task 
is to write down the equations for the singlet and triplet states separately. 
It is convenient to redefine the wave functions of the coupled states and the 
interaction term w by extracting the spin variables, that are included in the 
indices A, A' and so on: 

Xrj,s (k, p) = , 



^s,(T — s,a—s' ,s' 



(p, k - p; k - p', p') EE WA,K-A;K-A',A' 



Here as usual: K = {'k,a);A = {p,s);A = {p',s');a — 0, ±l;s — ±i. Lower 
(spin) indices of the interaction w are restricted to the values ±-i by definition. 

Let us rewrite Eq.([l^) using these notations for each spin projection cr of a 
coupled state. The equation for a = 1 becomes: 
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'-' 2 

+ i E "^hhhi (P, k - p; k - p', p') x^^ f (k, p') = 0. (16) 
p' 

This equation is apparently symmetrical in the spin indices and its solution 
should describe triplet states. We marked it by setting k = 1. 

The case ct = is more complex, because of the contributions from triplet 
and singlet states. Accordingly, one obtains the system of two equations: 

2Yl^wi _i._x x (p,k-p;k-p',p')4yHk,p') + 
p' 



wi _i.i _i (p,k-p;k-p',p')a;Q'^l"i(k,p')] = 0, 



[r!p + f^k-p - A„(7^)]a;^":"'(k, p) + 

i E[^-i,i;-ii (P,k - p;k - p', p') a;jT^(k,p') 
p' 

W_i 1.1 _i " "^^ " ^ ' 

2 ' 2 ' 2 ' 2 



(p, k - p; k - p', p') x[;'i (k, p')] = 0. (17) 



A nontrivial solution to this system exists only if both equations are equivalent. 
This equivalence is guaranteed by antisynimetrical properties of the x and w 
functions: 

4THk,p) = -4^:"^(k,k-p); 

Ws,s';s",s"' (P,P';P",P"') = -Ws',s;s",s"' (p', P; P", P'") = (P, P'; P'", P") 

For this reason we consider only the first equation of the system in the 
following. 

The equation for the symmetrical in spin indices solution (a;^^i"''(k, p') = 

u, 2 

^ri^!°i (k, P'))i the form: 
u, 2 

[f^p + f^k-p - UK)Wof (k, P) + 1 J^iii'i ,1 (p, k - p; k - p', p') + 

p' 

(p,k-p;k-p',p')]xjf)(k,p') =0 (18) 



Eqs.(16) and (18) describe the triplet states with different spin projection. 



It will be shown below that these equations are equivalent. 
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The equation for the antisymmetrical in spin indices solution to the system 
(ff^ (x[^°'i"^(k, p') = —x''i^'"l(k,p')) describes the singlet states: 

U, 2 ^1 2 

[f^p + f^k-p - Xo.iK)]xl^f^ (k, p) + i ^[wi 1 ^_ 1 , 1 (p, k - p; k - p', p') - 

p' 

(p,k-p;k-p',p')]xf,")(k,p') -0 (19) 



It is useful to express the interaction matrix w in Eqs.(p^,([T^) and ( |T9| ) in 
terms of the Hamiltonian parameters. In a more detailed form the Hamiltonian 
in the site representation may be written as: 



^ ^ ^ ] Uirin,mlnl^tn,s^ti,s'^rn' ,sCn' ,s (20) 
ranmfnf 

Tnnra'n' 



mnm'n' 



The first two terms are responsible for the triplet interaction (matrices C/^^ , uj^j^ 
are antisymmetrical) whereas the third term describes the singlet interaction 

(3) 

(matrix U„^^ „j,„/ is symmetrical). A comparison of the two equivalent expres- 
sions for the Hamiltonian (|l]), ( pO[ ) gives: 

If > , — U / / — r/^^^ 

m^s-.rLjS^m'^s-.njS — ^ m. — s\n.~s:m\ — s:n, — s — nin,nin' ^ 

TJ , , — TT I I ~ TI^^'' — r/^'^^ 

^ m,s:n, — s;m,s;n, — s — ^ m,s;n,s:m' , — s:n,s — ^ ran.mrnr ^ mn,mlnh 

Um, s-.n, -s-m' ,-s-n' ,s — Urn.-s-n.s-.rn' ,s:n' ,-s = C^mn,m'n' ~^ ^mn,m'n' ■ i"^^) 

Here the site and spin indices are shown separately for the sake of clarity. It 
seems reasonable to suppose that the form of the interaction depends only on 
the spin of a 4-fermion state, rather than on its spin projection. It leads to 
an equality C/^^^ = t^mn.m'n' ' that, being not very important, makes cal- 
culations more easy. Using Eqs.(^), (^2|) in the momentum representation we 
may insert the Hamiltonian parameters (i.e. matrixes U^^\ J7^^^ and C/*^'^-') into 
Eqs.(ll6l),(|l|) and (|l|) instead of the function w. After some transformations 



one obtains the equation for the singlet states: 

+ n^-p - \oA^)Wof\Kp) + (22) 



i(t/.^ -f t/.^)^c/S_^^,_^,,^,4"r)(k,p') = 0, . = ±1 
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Both of equivalent Eqs.(p^,(p^) for the triplet states may be transformed to: 
i^p + n^-p - Ai,„(k)]a;^;f ) (k, p) + (23) 

i(t/.^ +t/.^)5:c/«_^^,_^,,^,.(^f)(k,p') . 0, a = o,±i. 
p' 

Equations (^) , (^3|) may be solved numerically though they are too complicated 
yet. As it will be shown below, the account of the finite number of sites in 
the interaction parameters U'^^\ U^^\ C/^'^^ permits to simplify the calculations 
significantly. An example of such a model is discussed in the next section. 



2.3 Spectra of 2-particle excitations. 1 — 2 model 

Let us fix up the form of the fermion - fermion interaction. It takes into ac- 
count the interaction on the same or neighboring sites (1-2 model). The 
interaction parameters t/*^^', C/*^^', J/'-^-' (see expression (|20|)) in the momentum 
representation have the form: 



f^p!i-k;k-p'.p' = C4'i-k;k-p',p' = SwW'Ut2 EJcoKP - P', 1) - CO.s(k - P - p', 1)] 
C^S-k;k-p',p' - '^kk'{C^.l + Us2 E1[C0S(P - P', 1) + COs{k - p - p', 1)]}, 

where Ut2 is the interaction energy of the fcrmions in the triplet state and 
Us2,Usi are the intersite and onsite interaction of the fermions in the singlet 
state; 1 is a lattice vector. Substituting these expressions into Eqs.(^, [2^ ) one 
can obtain: 

x^if ) (k, p) - 2LkpUt2 J2 ■'^^"(^k - p, 1) ^ smi^k - p', l)xi]fHk, p') = (24) 
1 p' 

- for the triplet state and 



4°f ) (k, p) - Lfcp + 2(7.2 J2 - p, l)cos(ik - p', \)]xl^f (k, p') = (25) 

p' 1 

-for the singlet one. Here the notations are: 

i th{(3np/2) + th{(3n^_p/2) ^ 

Lkp t;—. 7r-\ — n o ' i^k = - M = f > cos(kl) - /z, 

^ '^(K,a)(.Kj - iip - iik-p J 
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p is the internal momentum variable and the hopping integral is supposed 
to be t(m — n) = t ^ Q only if m = n ± 1. These equations obviously lead to 
the following p - dependence of the pair wave function: 

dim 

x^^fHKp) =Lfcp^A(*)(k)sin(-k-p,l), (26) 

1=1 



dim 
1=1 



dim _j 

(k,p) =L,,[4^\k)+Y,4^\k)co8{-k-p,l)], (27) 



where index ^ = 1 for 1 = (;^, 0, 0), / = 2 for 1 = (0, ly,0), I S for 1 = (0, 0, 1,) 
and dim stands for the dimensionality of the system. The lattice vector 1 numer- 
ates in a natural way the triplet bands and therefore the number of the triplet 
bands is equal to the space dimensionality. The coefficients ^^^■'(k), ^[^■'(k) 
A^\k) are determined by substitution of these functional forms into Eqs.(p^, 
pq ) for the wave functions. The equation describing the triplet pairs appears to 
be diagonal after this substitution and is transformed into linear independent 
algebraic equations for the functions (k) : 

[1 - 2Ut2Hi{k)]A['\k) = 0,l< dim 

where Hi = J2p LkpSm'^i^k - p, 1). 

The excitation spectra are obtained from the equations: 

1 - 2Ut2Hi{k) = 0, (28) 
while the coefficients aI*-* (k) are derived from the normalization condition (^5|) 

which is now written as ^ xi^s"^*(k, p)xl^s" ^(k, p) = Sa.a'- With the ac- 
count of Eq.(|2^) one obtains: 

-1/2 



(^5:L?.,sin2(ik-p,l)j 



As a result of the diagonalisation the triplet wave functions ( pq ) appear to be 
rather simple: 

xiif)(k,p) =Lfe,pA[*)(k)sin(ik-p,l) (29) 

The singlet pair case is more complicated. Nevertheless one can obtain a linear 
algebraic system of 4 equations on the coefficients AQ^\k), ^|''-'(k) by substitut- 



ing expression (\21\ ) into Eq.(|25|): 

/ Ao \ 

Ai 
A2 

V ^3 ; 



(i - M) 



0, (30) 
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where M is the matrix of rank 4: 



M 



( Usllo Uslh Uslh Uslh \ 

2Us2h 2[/,2/ll 2Us2h2 2[/,2/l3 

2Us2h 2[/,2/21 2Us2h2 2f/,2/23 

V 2Us2h 2U,2hl 2Us2h2 2f/,2/33- / 



(31) 



and 



-^0 = Ep Lkp\ h = Ep cos(ik - p, 1); 
Iim = Ep ^fcp cos(ik - p, 1) • cos(ik - p, m). 

The matrix M is written for the 3 - dimensional system. The same matrix for 
the 2 - dimensional system may be obtained by omitting the last column and the 
last line. For a 1 - dimensional system this simple procedure must be repeated. 
The spectra of the singlet states are found by equating the determinant of the 
system (RQ) to zero: 



Det{i - M) = 0, 



(32) 



whereas the structure coefficients a\^^ (k) , (k) describing the wave function 
of the singlet pair ( p7| ) are found from Eq. ( |30| ) and the normalization condition: 



Now we have to find the excitation spectra of the local pairs by utilizing 
formulae ( psf ) - (^2|). The numerical calculations have been done in the units: 
Dpoi = 4. The results are represented in Fig.l for the interaction value of 
Usi = +8 (repulsive onsite interaction term); Us2 — — 6.8, J7t2 = —6.5 (attractive 
intersite terms for both the singlet and triplet states). 

Here the most interesting case is considered: onsite repulsion (Coulomb force 
predominance) and intersite attraction. The concrete value of onsite repulsion 
has only a small effect on the character of band structure. The explanation is 
clear: two coupled fermions are placed on neighbor sites mostly where attraction 
dominates and "don't feel" onsite repulsion. On the contrary, the value of the 
intersite attraction is important for the determination of low - lying bands. 
This situation has been analyzed in detail in The authors argued that the 
gap between the local pair band and the polaron band, that is placed above, 
increases with the coupling constant when considering onsite attractive Hubbard 
model. The substitution of the onsite attraction by the intersite one does not 
affect this qualitative result. Fig.l represents band structure of the singlet and 
triplet pairs at definite values of s-wave and p-wave attraction. The choice of 
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the attraction value Us2 = —6.8 and Ut2 = —6.5 gives reasonable mass of local 
pairs. 

It becomes evident that the band structure is not as simple as in the standard 
1-site attractive Hubbard model. There are additional gaps between the singlet 
and triplet bands, along with large ordinary gap. Each of the singlet and triplet 
band in its turn is splitted into two subbands . This could explain those contra- 
dictions that seem take place in determination of the gap value from different 
experiments. A one - particle (polaron) band is lying higher with the energy of 
its bottom being equal 0.97. It is necessary now to normalize our relative units. 
An appropriate value of a polaron bandwidth is about 50-lOOmeV. The undi- 
mentional polaron bandwidth is equal 4 in the units that were used throughout 
the paper. It means that our energy unit is about 20meV, The chosen value 
of the intersite attraction is 136mey(6.8 • 20meV), that seems reasonable. The 
value of the gap between the lowest local pair energy and lowest polaron energy 
depends strictly on the sort of pairing. Thus, in the case presented in Fg.l, the 
gapwidth between polaron band and s-pair band is about 207716^(0.97 • 20meV) 

Another interesting conclusion concerns a mass of the local pairs. According 
to the wide-spread point of view the local pair mass should be large enough 
this being considered as a shortcoming of local pair theories. Our calculations 
demonstrate (see Fig.l), that the situation is more favorable: the width of the 
lowest singlet band is only 3-4 times smaller, than the width of the fermion 
band. It gives the mass of the local pair m** ~ 3 47n*, where m* is the 
effective fermion mass. One can notice finally, that the triplet pairs appear to 
be several times heavier than the singlet ones for this set of parameters. 



3 Local pairs in external fields 



The excitation spectra and the local pair structure studied above may be used 
for theoretical investigation of thermodynamic and kinetic properties of the sys- 
tems with the strongly interacting electrons. On the other hand there are many 
experiments dealing with the interaction of HTSC materials with different ex- 
ternal fields. For the theoretical explanation of these experiments one has to 
know the form of the interaction of the pairs with the external fields. The 
corresponding vertex functions should be obviously dependent on the internal 
structure of the coupled state. The interaction of the fermion with the external 
field presented in Eq.(|l|) for the Hamiltonian is of rather general form. It in- 
cludes, for example, the interaction with neutrons and electromagnetic field. By 
integrating over the grassman variables in the same way as it was done before, 
one can obtain the action in the form (H), with the matrix Y depending on the 
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external field: 



{rpp' +ie5pp>)5^^' + B{e-e)pp> ij; ie + e 

Ypp' = \ / ,\ 

r[t + e j^^^ -{rpp, + ^e^ppO(5,,' - B{e - e ), 

In Eq.(|3^) all matrices are written in the momentum representation and an 
imaginary " time" Fourier - transformation is assumed. Expanding this expres- 
sion up to the second order in the fields ip, ip* and to the first order in the (weak) 
field B, one can obtain the effective action for the local pairs in the external 
field: 



Z = J DC*DCexp 5(C*, C), 
C) = d^"^* ik,E) {lE - A,„(k)) C^"") {k,E) 

{k,E) rt?''^'"') (k,i?, k',E') cir'"') {k',E') . 

Here the summation on the repeated indices is assumed. Eq.(p^ and the or- 
thonormal basis (|l3|) were used in the derivation of this equation. The vertex F 
has the form: 

fI:?-'"') (k,E,k',E') 25:i?.,_,._.(k-k',i?-ij')Ftr'"')(k-k'). 

S 

Here the formfactor 

^i:r'"'^(k-kVE-^"-:(k.k-p)4:;^^^ (35) 
p 

depends only on the structure of the coupled state. The external field is written 
in the more convenient form - with the spin index extracted and translational in- 
variance taken into account: Bpp> (E) = Bg^.s^, (p — p', E). As it was mentioned 
above, the vertex of the interaction of pairs with the external field depends on 
the frequency, the wave vector of the external field and the wave functions of 
the fermion pairs. One can see that the slowly varying field is not sensitive to 
the pair structure. A simple test shows that in this case the vertex of the singlet 
local pair interaction with the external field equals to zero. For the triplet pairs 
it is twice the interaction of the fermion with the field. The behavior of formfac- 



tor (35) plays a significant role in the theoretical explanation of experiments on 
the scattering of particles on the local pairs. For example, at low temperatures 
and small transferred momenta the transitions appear to be suppressed. The 
reason is quite obvious: at the low temperature the local pairs being bosons are 
condensed at the lowest level with the momenta k = 0. The transitions into the 
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upper bands with approximately the same momentum are suppressed due to the 
orthogonality of the wave functions with different energies: F(k — k' = 0) = 0. 

The typical formfactors of the transitions from the lowest singlet band to the 
upper singlet or triplet one are presented in Fig. 2. There is also the formfactor 
describing the transitions inside the lowest singlet band. The transitions into the 
second triplet band arc forbidden for chosen transferred momenta k = {k^^O) 
due to the orthogonality of the wave functions. The evident momentum de- 
pendence of the formfactor confirms its importance in the interpretation of 
experiments. 

4 Conclusion 

The mathematics of the local pair theory was developed. The ad hoc calculation 
of band structure of the local pairs and their interaction with external fields are 
presented here. It is shown that numerical calculations based on analytical 
formulas and being relatively simple lead to exact qualitative predictions. The 
method can be easily generalized. The momentum dependency of the vertex is 
essential for theoretical explanation of experiments . 
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